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ABSTRACT 

We examine Lyman continuum extinction (LCE) in H II regions by comparing infrared fluxes of 49 H 
II regions in the Galaxy M31, M33, and the LMC with estimated production rates of Lyman continuum 
photons. A typical fraction of Lyman continuum photons that contribute to hydrogen ionization in the 
H II regions of three spiral galaxies is < 50 %. The fraction may become smaller as the metallicity (or 
dust-to-gas ratio) increases. We examine the LCE effect on estimated star formation rates (SFR) of 
galaxies. The correction factor for the Galactic dust-to-gas ratio is 2-5. 

Subject headings: dust, extinction — H n regions — galaxies: ISM — infrared: ISM: continuum — 
stars: formation 



1. INTRODUCTION 

Dust grains exist everywhere — from circumstellar en- 
vironment to interstellar space. Hence, radiation from ce- 
lestial objects is always absorbed and scattered by dust 
grains. Without correction for the dust extinction, we in- 
evitably underestimate the intrinsic intensity of the radia- 
tion, and our understanding of physics of these objects 
may be misled. Thus, dust extinction for nonionizing 
photons (A > 912 A) in the interstellar medium (ISM) 
has been well studied to date (e.g., interstellar extinction 
curve: Savage & Mathis 1979; Seaton 1979; Calzetti, Kin- 
ney, & Storchi-Bergmann 1994; Gordon et al. 2000). 

In H II regions, Lyman continuum (LC) photons also 
suffer extinction by dust grains (e.g., Ishida & Kawajiri 
1968; Harper & Low 1971). We refer to this another type 
of extinction as the Lyman continuum extinction (LCE). 
Indeed, the fraction of LC photons absorbed by dust in 
H II regions can be large (e.g., Petrosian, Silk, & Field 
1972; Panagia 1974; Mezger, Smith, & Churchwell 1974; 
Natta & Panagia 1976; Sarazin 1977; Smith, Biermann, & 
Mezger 1978; Aannestad 1989; Shields & Kennicutt 1995; 
Bottorff et al. 1998). Petrosian et al. (1972), for example, 
have estimated the fraction of LC photons contributing to 
hydrogen ionization to be 0.26 in the Orion nebula. 

If a large amount of LC photons is absorbed by dust in 
H II regions, the effect of the LCE on estimating the star 
formation rate (SFR) will be very large. This is because we 
can only estimate the number of ionizing photons, LC pho- 
tons used for hydrogen ionization, from the observation of 
hydrogen recombination lines or thermal radio continuum. 
When we define the fraction of ionizing photons estimated 
from observations as /, we obtain 

Kc = fN LC , (1) 
where -ZVlc and iV LC are the intrinsic and apparent produc- 
tion rates of LC photons, respectively. Unless we correct 
the observed data for the LCE, we cannot obtain the ac- 
tual LC photon production rate, and then, the true SFR. 
Therefore, we should estimate the fraction, /, in H II re- 
gions. However, the LCE has not been discussed in the 
context of estimating the SFR of galaxies before Inoue, 

See also Inoue, Hirashita, & Kamaya (2001b) 



Hirashita, & Kamaya (2001a, hereafter Paper I). 

Paper I has proposed two independent methods for es- 
timating /, and applied the methods to the Galactic H n 
regions. Then, they have shown that there are a number 
of H II regions with / < 0.5 in the Galaxy. Moreover, 
they have found that / can be approximated to be only a 
function of dust-to-gas ratio. If the approximation can be 
applied to the nearby spiral galaxies and their observed 
dust-to-gas ratios are adopted, / averaged over galactic 
scale is about 0.3. Therefore, the SFR corrected by the 
LCE effect is 3 times larger than the uncorrected SFR. 

The aim of this paper is that we confirm these results 
of Paper I not only for H n regions in our Galaxy but for 
extragalactic H n regions. We also discuss if the effect of 
the LCE on estimating the SFR of galaxies is really impor- 
tant. We describe the method for estimating / in section 
2. Then, / is determined for the H II regions in the local 
group galaxies in section 3. We discuss the dependence of 
/ on the dust-to-gas ratio in section 4, and our conclusions 
are summarized in the final section. 

2. METHOD FOR ESTIMATING THE EFFECT OF LCE 

To estimate / in equation (1) for each H n region, the 
apparent production rate of LC photons in each region, 
Y LC , is compared to its IR luminosity, Lir. In this paper, 
we use a method for estimating / presented in section 3 of 
Paper I. 

According to equation (12) in Paper I, we obtain 



J 0.28 + 5.6(^7/^,60) ' 

where e denotes the average efficiency of dust absorption 
for nonionizing (A > 912 A) photons, and Lir. 7 and -/V LC 50 
are the IR luminosity and LC photon production rate nor- 
malized by 10 7 Lq and 10 50 s _1 , respectively. The equa- 
tion is based on the theory of the IR emission from H 11 
regions by Petrosian et al. (1972) and the following dis- 
cussion by Inoue, Hirashita, & Kamaya (2000) : . It is also 
worthwhile to note that Lir^ /N^ c 50 is proportional to the 
flux ratio of the IR emission to the thermal radio emission, 
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Ha line, etc. If we express in the form of the flux ratio, 
the determination of / is free from the uncertainty of the 
distance to the object. 

Here, wc describe some assumptions adopted in the 
derivation of equation (2). Since Lir is the total energy 
emitted by dust, the wavelength range is set to 8-1000 fim, 
which covers almost all wavelength range of dust emission. 
Also, it is supposed that the IR radiation from H II regions 
is dominated by large (~ 0.1/xm) grains in thermal equi- 
librium with ambient radiation field. In other words, we 
neglect the contribution to the total IR luminosity of very 
small grains (< O.Ol^im) and polycyclic aromatic hydro- 
carbons (PAHs), which emit mainly in the mid-infrared 
(MIR, 8-40 fjaa) (e.g., Dwek et al. 1997). In principle, 
the assumption would make us underestimate the total 
IR luminosity of dust in H II regions because such MIR 
emission is observed from these regions (e.g., Pi§mi§ & 
Mampaso 1991; DeGioia-Eastwood 1992; Leisawitz et al. 
1998). Fortunately, it seems that the contribution of the 
MIR luminosity to that of the total IR is enough small 
to be neglected (Figure 1 in Leisawitz et al. 1998). Thus, 
the IR spectral energy distribution (SED) of H n regions 
is assumed to be a modified black-body function with one 
temperature through the paper. 

Moreover, we do not take account of the existence of he- 
lium. This is because almost all (96% according to Mathis 
1971) photons produced by helium recombination can ion- 
ize neutral hydrogen. We do not consider the escape of LC 
photons from H n regions, either. Indeed, less than 3 % 
of LC photons can escape from the star-forming regions in 
nearby galaxies (Leitherer et al. 1995). In short, we deal 
with an ideal H n region that consists of only hydrogen 
and dust, and where no LC photons leak. 

Let us determine e, which denotes the average efficiency 
of dust absorption for nonionizing photons. The total (ab- 
sorption and scattering) optical depth of dust, t\, is ex- 
pressed by 

abs A . 

1 - lu\ 2.5 loge 

where r^ bs , ui\, and A\ are the absorption optical depth 
of dust, the dust albedo, and the amount of extinction at 
wavelength, A, respectively. We can determine A\ from 
an extinction curve, X\, and its normalization, Eb-v, i-e. 
X\ = A\/E B _y The efficiency of dust absorption at 
wavelength, A, e\ is defined as 

e x = 1 - e- r " s = 1 - 10 -oax x e B - V (i-u x ) (4) 

Then, averaging e\ over A with weight of the intrinsic stel- 
lar flux density yields the average efficiency, e. 

We show e as a function of Eb-v in Figure 1. In this 
calculation, we set the stellar spectrum to be simply the 
black-body of 30000 K. The change of the temperature 
dose not alter the result significantly. For example, if we 
choose 50000 K as the effective temperature of stars, the 
value of e increases with a factor of about 1 %. The calcu- 
lation of the average is performed from 1100 to 3500 A. We 
adopt two types of the interstellar extinction curve, X\: 
one is that in our Galaxy (Seaton 1979), the other is that 
in the Large Magellanic Cloud (LMC) (Howarth 1983). 
Each case is shown in Figure 1 by the solid and dashed 

2 The estimated production rate of LC photons of M8, Orion, and 
references are different from those used in this study. 



lines, respectively. Moreover, we set the dust albedo to 
be a constant of 0.5 for simplicity, although it varies with 
wavelength (Witt & Lillie 1973; Lillie & Witt 1976). The 
results in this case are shown by the thick lines. Also, the 
results without scattering (i.e. u\ = 0) are shown by the 
thin lines for comparison. 

Once we determine e for H n regions from Figure 1, we 
can estimate / from the ratio of the IR luminosity to the 
observed LC photon production rate via equation (2). In 
the next section, we will examine / for H n regions in some 
Local Group galaxies. 

3. FRACTION / OF H II REGIONS IN THE LOCAL GROUP 
GALAXIES 

In order to determine /, the data of the IR luminos- 
ity and the production rate of LC photons for the indi- 
vidual H II region are required. But, when we compare 
two photometric data, we must match their apertures and 
resolutions with each other, because, in general, different 
observations have different aperture and resolution. It pre- 
vents us from constructing uniform data of a large number 
of H II regions. These effects produce one of the largest 
uncertainties in the current analysis. 

We have gathered the data that allow us to estimate 
/ with proper accuracy, though the sample size is rather 
small. In this section, we estimate / of H n regions in our 
Galaxy, M31, M33, and LMC. 

3.1. Our Galaxy 

Suitable data for our analysis have been compiled by 
Wynn- Williams & Becklin (1974). The data are read out 
from their figure 9 directly, and are tabulated in Table 
1. The column 1 is the name of H n regions. Since two 
separate components of W3, W51, and NGC 6537 are plot- 
ted in figure 9 of Wynn- Williams & Becklin (1974), wc list 
both components of these regions in Table 1. The columns 
2 and 3 are the logarithmic production rate of LC photons 
estimated from radio observations 2 and FIR (40-350 /im) 
luminosity, respectively. The FIR (40-350 /im) luminosity 
is almost equal to the total IR (8-1000 /im) luminosity if 
the dust SED is assumed to be the modified black-body 
function with the temperature of 30 K and the emissivity 
index of 1. 

Since Eb-v ~ 1 for the Galactic H n regions (e.g., Ca- 
plan et al. 2000), Figure 1 tells us that e ~ 1, even in the 
case with the scattering for nonionizing photons. Thus, 
we can safely assume e = 1 for the Galactic H n regions, 
so that we can determine / of these regions via equation 
(2) (column 5). The mean value of / is 0.45. Thus, we 
find that about half of LC photons is absorbed by dust in 
H II regions. This is one of the results in Paper I. 

3.2. M31 

Xu & Hclou (1996) have studied the properties of dis- 
crete FIR sources in M31 by using the IRAS HiRes image 
(Aumann, Fowler, & Melnyk 1990; Rice 1993). Thirty nine 
sources are extracted from the 60 fim map. They have 
also examined the correlation between these FIR sources 
and Ha sources presented by Walterbos & Braun (1992). 

in this paper differ from those in Table 1 of Paper I because their 
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Since Waltcrbos & Braun (1992) have surveyed only the 
northeast half of M31, 12 out of 39 FIR sources have Ha 
counterparts. These twelve FIR sources make the sample 
for the current analysis. In Table 2, some properties of the 
sample H n regions in M31 are listed. 

The column 2 of Table 2 is the measured flux density at 
60 /mi. Unfortunately, only the 60 /xm flux densities of al- 
most all sample regions are available. Hence, we estimate 
total IR luminosity of dust emission in each region from the 
60 /im flux density by the following way: Assuming the IR 
SED of these FIR sources to be the modified black-body 
function with temperature of 30 K and emissivity index 
of 1, we obtain L m /L Q = 3.50 x (D/kpc) 2 (F 60 /mi/Jy), 
where the wavelength range of the IR luminosity, £ir, is 
8-1000 /im, and D = 760 kpc is the adopted distance to 
M31. The derived IR luminosities are shown in the col- 
umn 3 of Table 2. The adopted dust temperature of 30 
K is reasonable because dust temperature of 10 isolated H 
II regions whose flux densities at both 60 and 100 /<m are 
listed in Table 3 of Xu & Hclou (1996) is estimated to be 
29 K on average. 

Since the angular resolution of Ha image is much higher 
than that of 60 image, numerous Ha sources are in- 
cluded in the area over which the flux density of one FIR 
source at 60 //m is integrated. The sum of the Ha fluxes 
in the area of each FIR source is shown in the column 4 
of Table 2. These values are corrected for the [Nil] con- 
tamination by assuming a constant ratio of [Nil] /Ha = 0.3 
(Walterbos & Braun 1992). 

When we estimate the production rate of LC photons 
from Ha flux, we must correct it for the interstellar extinc- 
tion. Walterbos & Braun (1992) have commented that the 
mode of the interstellar extinction at Ha is 0.8 mag. Since 
this value of extinction is estimated from the column den- 
sity of Hi gas, it is a typical extinction averaged over the 
ISM of M31. If we assume the Case B and the gas temper- 
ature of 10000 K (Osterbrock 1989), we obtain iV^/s^ 1 = 
8.77x 10 55 (Z?/kpc) 2 (F HQ /ergss- 1 cm" 2 ) x 10°- 4Ah «, where 
Fh q is Ha flux (column 4) and A^ a = 0.8 mag is the as- 
sumed amount of extinction at Ha line. The derived N^ c 
is in the column 5. The ratio of the IR luminosity to the 
LC photon production rate is shown in the column 6. 

Now, we need to estimate e. Adopting Ay /Eb-v = 3.1 
and the Galactic extinction curve, Eb-v is ~ 0.3 mag 
when An a ~ 0.8 mag. Then, we can assume e = 0.7 for 
the H II regions in M31 from the line for the dust albedo 
of 0.5 in Figure 1. Finally, we determine / via equation 
(2) (column 7). The mean / is 0.38, which is almost equal 
to that of the Galactic H II regions. 

Also, Devereux et al. (1994) have compared the Ha 
emission map of M31 with the IR emission map produced 
by the maximum correlation method (IRAS HiRes im- 
age). They divided the Ha emission into three compo- 
nents; the star-forming ring, the bright nucleus, and diffuse 
(unidentified) component. The star- forming ring, whose 
diameter is 1?65, has the observed Ha+[Nll] luminosity 
of 4.77 x 1O 6 L . They estimated its IR (8-1000 /mi) lu- 
minosity to be 1.08 x 1O 9 L0 by assuming the modified 
black-body function with 26 K and the index of 1. If we 
regard the star-forming ring as an aggregation of a lot of 
H II regions, and assume A^ a = 0.8 mag and [Nn]/Ha = 
0.3, we obtain ^ir, 7/^0,50 ~ 0-5 for the ring. This cor- 



responds to / ~ 0.3 when e = 0.7. The value of / shows 
a good agreement with the mean / of the sample in Table 
2. 

3.3. M33 

We find suitable data of H 11 regions in M33 in Devereux, 
Duric, & Scowen (1997), who examine correlation between 
Ha, FIR, and thermal radio emission of M33. They have 
reported some properties of 8 isolated H II regions. We 
determine / of these regions. The sample properties are 
tabulated in Table 3. 

In the column 2, FIR (40-1000 /tm) luminosities of the 
H II regions, which are determined from the flux densi- 
ties at 60 and 100 /mi measured from IRAS HiRes images, 
are presented. In the column 3, the 5 GHz flux densi- 
ties (see also Duric et al. 1993) are shown, and the esti- 
mated LC photon production rates are in the column 4 
via Nlc/s- 1 = 8.88 x 10 46 ( J D/kpc) 2 (F 5GH z/Jy) (Condon 
1992). For the distance to M33, 840 kpc is adopted. If 
dust SED can be fitted by modified black-body function 
(emissivity index of 1, 30K), we have I/fir (40-1000 /an) 
~ Lir (8-1000 /mi). We determine the ratio of the IR 
luminosity to the photon production rate under this as- 
sumption (column 5). 

Devereux et al. (1997) have also commented on the prop- 
erties of the interstellar extinction for the H 11 regions in 
M33. They have compared Ha emission with thermal ra- 
dio emission, and argued in their figure 7 that the average 
extinction is Ay <~ 1 mag (i.e. iia ~ 0.8 mag), although 
there is a large dispersion corresponding to ±1 mag. Thus, 
we can assume e = 0.7 for all the sample regions in M33 
as well as those in M31. The determined values of / are 
listed in the column 6 of Table 3. The mean is 0.41, al- 
though the large dispersion of the interstellar extinction 
may cause a large dispersion of /. 

Thus, we saw that / ~ 0.4 or less for the H II regions 
in M31 and M33. Hence, we conclude that the typical 
values of / in these galaxies are nearly equal to that in 
our Galaxy. Therefore, we suggest that the fraction of LC 
photons consumed by dust in the H 11 regions of other spi- 
ral galaxies may be a > 50% as well as that in the Galaxy, 
M31, and M33. 

3.4. LMC 

DcGioia-Eastwood (1992) have provided the data of 6 H 
11 regions in the LMC. The IR flux densities of the sample 
regions have been measured from the co-added survey data 
by IRAS. Relatively isolated regions have been chosen as 
the sample regions in order to perform a reasonable sub- 
traction of the background flux. DeGioia-Eastwood (1992) 
have reported the IR, Ha, and radio fluxes of the sam- 
ple regions whose apertures have been matched carefully 
among images of different wavelengths. The properties of 
the data are shown in Table 4. 

The columns 2 and 3 are IRAS 60 and 100 /im flux 
densities, respectively. The IR (8-1000 /mi) luminos- 
ity is calculated from these flux densities via Lir/L q = 
0.61 x (D/kpc) 2 [2.58(F 60 /Jy) + (F 100 /Jy)]. Here, we as- 
sume the dust SED to be modified black-body function 
with 30 K and index of 1, that is, L m (8-1000 /xm) ~ 1.66 
-Lfir (40-120 /xm). Then, the estimated total IR lumi- 
nosities are shown in the column 4. The adopted distance 
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to the LMC is 51.8 kpc. From the flux densities at 5 GHz 
in the column 5, the production rates of LC photons are 
calculated via Condon's formula in section 3.3 (column 6). 

The color excess of each region, Eb-v, estimated from 
the observed Balmer decrement and corrected for a fore- 
ground reddening is listed in the column 8. Using the line 
of the LMC extinction law and the dust albedo of 0.5 in 
Figure 1, we obtained the value of e from each Eb-v (col- 
umn 9). Finally, / is determined for each sample region 
(column 10). The mean / is 0.74. This is obviously larger 
than those of the Galaxy, M31, and M33. 

DeGioia-Eastwood (1992) have also determined / for 
these regions. Her mean / is 0.65 ± 0.13, which is consis- 
tent with that in the current paper. Here, we should note 
that / of DeGioia-Eastwood (1992) is determined from the 
dust optical depth for UV (1000-3000 A) photons. By the 
definition of /, however, it should be determined from the 
dust optical depth for LC photons (A < 912 A) such as we 
stated in Paper I or this paper. 

On the other hand, Xu et al. (1992) have compared the 
IR, Ha, and radio emission from the LMC, and have ex- 
amined more global correlations among these emissions. 
They have displayed appealing contour maps of these 
wavelengths, which strongly make us believe that these 
three kinds of emission have the same origin, i.e. the star- 
forming regions. Indeed, there are excellent coincidences 
among IR, thermal radio, and Ha sources (see also Mai- 
hara, Oda, & Okuda 1979; Fiirst, Reich, and Sofue 1987 
for the Galactic sources). According to Xu et al. (1992), 
the FIR (40-120 /zm) flux associating with the star forma- 
tion in the LMC is 28.9 x 10~ 10 W mT 2 , and the thermal 
radio flux density at 6.3 cm in the same area is 131 Jy. 3 
Since the corresponding Li^/N^ c is 0.13, / ~ 0.7 when 
we assume e = 0.4 (the average value for the sample in 
Table 4). This is in good agreement with the mean / of 
the sample in Table 4. 

Therefore, we conclude that / of the H n regions in 
the LMC is about 0.7 and larger than those in the spiral 
galaxies like the Galaxy, M31, and M33. The larger / may 
be caused by the lower metallicity of the LMC. A lower 
mctallicity is likely to lead to a smaller dust-to-gas ratio, 
so that the fraction of ionizing photons, /, increases. The 
issue will be discussed more closely in the next section. 

4. DISCUSSIONS 

We have estimated / for each H n region in some Local 
Group galaxies by comparing its observed IR luminosity 
with production rate of LC photons. Here, we discuss what 
the determined / values suggest. 

4.1. Value of f as a function of metallicity /dust-to- gas 

ratio 

First, we discuss the relationship between the mean / 
of each galaxy and its metallicity (or dust-to-gas ratio). 
We naturally expect that the effect of the LCE becomes 
larger as the dust content increases. Then, we also expect 
the mean / of sample H II regions to be a function of dust- 
to-gas ratio or metallicity of their host galaxy. In order to 



examine the issue, we summarize some parameters of each 
galaxy in Table 5. 

The mean / with its sample standard deviation is shown 
in the column 2 of Table 5. Also, we find metallicities and 
dust-to-gas mass ratios of sample galaxies in the columns 
3 and 4, respectively. Here, the metallicities arc taken 
from van den Bergh (2000). The dust-to-gas mass ratios 
except for that of M31 are taken from Issa, MacLaren, 
& Wolfcndale (1990), who have estimated the dust-to-gas 
ratio from the ratio of the amount of visual extinction to 
the surface density of H I gas. For M31, we determined 
the dust-to-gas ratio from the observed dust mass and gas 
mass directly. According to the observation of Haas et 
al. (1998) by ISO, M31 have 3 ± 1 x 10 7 M Q as the dust 
mass. Also, M31 contains H I gas of 3.8 x 10 9 M© (Cram, 
Roberts, & Whitehurst 1980) and H 2 gas of 2.7 x 1O 8 M 
(Dame et al. 1993). Thus, the dust-to-gas mass ratio of 
M31 is 7.4 x 10~ 3 . Moreover, the dust-to-gas ratios of sam- 
ple galaxies are normalized by a typical one of the ISM in 
our Galaxy, 6 x 10~ 3 (Spitzer 1978). 

In Figure 2, mean / values of H II regions in each galaxy 
are plotted as a function of metallicity. We may find a 
trend that the mean / decreases as the metallicity in- 
creases. We also show that mean / as a function of dust- 
to-gas mass ratio in Figure 3. In this figure, the mean / 
seems to become small as the dust-to-gas ratio becomes 
large. However, since the number of sample galaxies is 
still small, more investigation is needed. For example, the 
dependence of / on mctallicity will become clearer if we 
examine H II regions in the Small Magellanic Cloud or 
other metal poor galaxies, or examine the relation of this 
issue to the abundance gradient of galaxies. 

Let us see how to produce two solid lines in Figure 3. 
In Parer I, we have formulated the dust optical depth for 
LC photons, Td, as a function of its dust-to-gas ratio (see 
also Hirashita et al. 2001); 

Td = x \rT~l ' ( 5 ) 

where Td is evaluated over the actual ionized radius of 
a spherical H n region containing dust grains uniformly. 
Here, we estimate Td at the Lyman limit approximately. 
In the right hand side, V is the dust-to-gas mass ra- 
tio which is normalized by a typical Galactic ISM value, 
Pmw = 6 x 10~ 3 (Spitzer 1978). The coefficient, x, is 
a factor depending on the production rate of LC photons 
and the gas density of the H II region. For a galaxy having 
the Galactic dust-to-gas ratio, the coefficient, x, means the 
dust optical depth for LC photons itself. We find, in Paper 
I, x ~ 2 for the observed LC photon production rates and 
electron number densities in the sample of Galactic H n 
regions. 4 Petrosian et al. (1972) have derived the relation 
between the monochromatic dust optical depth, Td, and / 
as the following; 



/ = 



(6) 



3{e-d( T 2 -2T d + 2) -2} • 
The relation is shown in figure 1 of Paper I. When T> is 
given, therefore, we determine / via equations (5) and (6). 
The data points in Figure 3 are well reproduced by setting 
x ~ 1-2 in equation (5). 

3 Their integrated area is not whole of the LMC. See their description for detail. 

4 This sample of Galactic H II regions differs from the sample in Table 1 of this paper, and the method estimating / is also different from that 
in the paper. See Paper I for detail. 
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Here, wc assume that the dust-to-gas ratio in H n re- 
gions to be almost equal to a typical value of that in the 
ISM of the host galaxy, although the ratio in local may 
vary significantly across the galaxy (Stanimirovic et al. 
2000). On the other hand, the mean value of the dust-to- 
gas ratio of the H n gas in our Galaxy is almost equal to 
those of the H I and H 2 gas (Sodroski et al. 1997). Thus, 
the assumption may be valid globally. 

In this paper, the determined / values are somewhat 
larger than those in Paper I systematically. That is, the 
suitable x parameter in equation (5) for the data points is 
1 rather than 2 which is the best fit value in Paper I. This 
may be because we neglect the MIR component emitted 
by very small grains and PAHs in the dust IR luminosity 
(see section 2). Thus, / determined here may be an upper 
limit. Also, / in Paper I is considered to be a lower limit 
because the scattering of LC photons by dust is not taken 
into account in Paper I. 

Anyway, the model presented by equations (5) and (6) 
can reproduce the trend of the observational data, when 
the coefficient, x, is set to between 1 and 2. In short, the 
V-f relation falls in the area between two curves in Figure 
3. Therefore, we conclude that / is expected a function 
of dust-to-gas ratio (or metallicity) : as the dust content 
increases, / becomes smaller. 

4.2. Effect of LCE on estimating SFR 

Finally, we discuss the effect of the LCE on estimating 
the SFR. By the definition of / in equation (1), wc obtain 
the correction factor for the SFR by the LCE as 1//. If 
/ is a function of dust-to-gas ratio, the correction factor 
is also a function of dust-to-gas ratio. This is shown in 
Figure 4. 

Clearly, we find that the correction factor increases 
as dust-to-gas ratio increases. For nearby spiral galax- 
ies, their dust-to-gas ratios distribute around the Galactic 
value (Alton et al. 1998; Stickel et al. 2000, see also Pa- 
per I). Indeed, even the starburst galaxies (e.g., M82) and 
the ultra- luminous IR (ULIR) galaxies (e.g., Arp 220) are 
likely to have about the Galactic dust-to-gas ratio (Kriigcl, 
Steppe, & Chini 1990; Lisenfeld, Isaak, & Hills 2000). For 
such dust-to-gas ratio, we expect that the correction fac- 
tor is 2-5. Thus, the effect of the LCE is more important 
than the uncertainty of the IMF, which is about a factor 



of 2 (e.g., Inouc et al. 2000). Therefore, we should take 
account of the LCE effect when we estimate the SFR of 
nearby spiral galaxies. 

5. CONCLUSIONS 

To examine the Lyman continuum extinction (LCE) in 
H II regions, we compared the observed infrared emission 
with production rate of Lyman continuum photons for 49 
H II regions in our Galaxy, M31, M33, and LMC. Then, 
we estimate the fraction of Lyman continuum photons con- 
tributing to hydrogen ionization, /, in these regions. We 
reached the following conclusions: 

[1] In many H II regions, / is smaller than 0.5. The mean 
/ of the sample regions in spiral galaxies (our Galaxy, M31, 
and M33) is about 0.4. On the other hand, the mean / in 
the H II regions in the LMC is about 0.7. 

[2] The mean / of sample H II regions in each galaxy 
may be a function of metallicity or dust-to-gas ratio of 
their host galaxy. Namely, / decreases as the metallicity 
or dust-to-gas ratio increases. 

[3] The observational trend is reproduced very well by 
the model that the dust optical depth for Lyman contin- 
uum photons over the ionized region is proportional to the 
dust-to-gas ratio. The dust optical depth is in between 1 
and 2 for the Galactic dust-to-gas ratio. 

[4] It is expected that the correction factor for the star 
formation rate by the LCE increases as dust-to-gas ra- 
tio (metallicity) increases. The expected correction factor 
for the galaxies with the Galactic dust-to-gas ratio is 2-5. 
Therefore, the LCE effect should be taken into account 
when we estimate the SFR of the nearby spiral galaxies, 
since the dust-to-gas ratios of these galaxies are almost 
same as that in the Galaxy. 
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Table 1 
The Galactic H ii regions 





W N' 


log i^FIR 


L ffi/ JV LC 






l b J 


IT \ 


±JQ j ±U b ) 




n ^ 
UJ 


(91 


\ 6 ) 


(A\ 

I 4 ; 




NGC2024 


47.7 


4.3 


0.40 


0.40 


W3 


48.2 


5.3 


1.26 


0.14 




49.4 


6.1 


0.50 


0.32 


M8 


48.3 


4.7 


0.25 


0.59 


NGC6357 


48.4 


5.2 


0.63 


0.26 




48.4 


5.0 


0.35 


0.44 


IC4628 


48.6 


5.4 


0.63 


0.26 


Orion 


48.7 


5.2 


0.35 


0.44 


G343.4-0.4 


48.7 


5.5 


0.56 


0.29 


DR15 


48.8 


5.8 


1.12 


0.15 


G5.9-0.4 


49.1 


5.6 


0.32 


0.49 


W75/DR21 


49.2 


5.3 


0.13 


1.01 


RCW117 


49.5 


6.1 


0.35 


0.44 


G351. 6+0.2 


49.7 


6.1 


0.25 


0.59 


W58 


49.7 


6.0 


0.20 


0.71 


G351. 6-1.3 


49.8 


6.3 


0.32 


0.49 


W51 


49.8 


6.2 


0.25 


0.59 




50.4 


6.9 


0.32 


0.49 


RCW122 


49.9 


6.3 


0.25 


0.59 


SgrC 


50.0 


6.8 


0.63 


0.26 


M17 


50.1 


6.5 


0.25 


0.59 


Sgr B 


50.5 


7.1 


0.35 


0.44 


W49 


50.7 


7.3 


0.35 


0.44 



Note. — Col. (1): Object name. Two separate components 
in W3, W51, and NGC6357 are both listed. Col. (2): LC photon 
production rates estimated from radio observations. Col. (3): Ob- 
served FIR luminosities. Col. (4): Ratios of IR luminosity to LC 
photon production rates. Col. (5): Estimated / from equation (2) 
by assuming e = 1. 
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Table 2 



H II REGIONS IN M31 



Object ID 










W^LC 


/ 




(Jy) 


(1O 6 L ) 


(10~ 12 erg s _1 cm~ 2 ) 


(io 5 ^- 1 ) 


(IO^q/IO 5 ^- 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


16 


3.50 


7.07 


3.33 


3.53 


0.20 


0.59 


17 


3.40 


6.87 


1.81 


1.92 


0.36 


0.36 


22 


4.56 


9.21 


4.37 


4.63 


0.20 


0.59 


23 


4.53 


9.15 


2.23 


2.36 


0.39 


0.34 


25 


3.96 


8.00 


0.77 


0.82 


0.98 


0.14 


27 


2.58 


5.21 


1.83 


1.94 


0.27 


0.46 


28 


6.55 


13.2 


3.51 


3.72 


0.36 


0.36 


30 


1.84 


3.72 


0.69 


0.73 


0.51 


0.26 


31 


6.52 


13.2 


2.08 


2.20 


0.60 


0.23 


33 


7.12 


14.4 


3.58 


3.79 


0.38 


0.34 


35 


2.27 


4.59 


1.12 


1.19 


0.39 


0.34 


39 


2.81 


5.68 


2.33 


2.47 


0.23 


0.53 



Note. — Col. (1): Identification of 60 sources in Xu & Hclou (1996). Col. (2): Extracted 60 
/mi fluxes from IRAS HiRes image by Xu & Hclou (1996). Col. (3): Estimated IR luminosities by 
adopting the distance of 760 kpc. Col. (4): Ha fluxes corresponding with 60 fim sources. Col. (5): LC 
photon production rates estimated from Ha fluxes in col. (4) by adopting the distance of 760 kpc and 
Ana = 0.8 mag. Col. (6): Ratios of IR luminosity to LC photon production rates. Col. (7): Estimated 
/ from equation (2) when we assume e = 0.7. 
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Table 3 



H II regions in M33 



Object 


Lfir 


-F5GHZ 




W^LC 

(IO^q/IO 5 ^- 1 ) 


f 




(10 7 L Q ) 


(mJy) 


(KPs- 1 ) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


VGHC 20,25 


0.56 


0.8 


0.50 


1.12 


0.13 


IC 131 


1.00 


0.7 


0.44 


2.27 


0.064 


IC 133 


1.67 


6.6 


4.14 


0.40 


0.33 


NGC 595 


0.79 


18.0 


11.3 


0.070 


1.23 


VGHC 46,52 


1.34 


1.8 


1.13 


1.19 


0.12 


VGHC 104 


0.58 


1.1 


0.69 


0.84 


0.17 


VGHC 97,98 


0.47 


4.5 


2.82 


0.17 


0.67 


NGC 604 


7.50 


60.0 


37.6 


0.20 


0.59 



Note. — Col.(l): Object name. Col. (2): FIR luminosities determined from 
IRAS HiRes images of 60 and 100 /im by adopting the distance of 840 kpc. 
Col. (3): Observed 5 GHz flux densities. Col. (4): LC photon production rates 
estimated from col. (3) by adopting the distance of 840 kpc. Col. (5): Ratios 
of IR luminosity to LC photon production rates. Col. (6): Estimated / from 
equation (2) by assuming e = 0.7. 
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Table 4 
H n regions in LMC 



Object 






Lm 


^5GHz 


N Lc 


W^LC 


Eb-v 


e 


/ 




(Jy) 


(Jy) 


(1O 6 L ) 


(Jy) 


(io 5 ^- 1 ) 


(lO^o/lO 5 ^" 1 ) 


(mag) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


MC 18 


1229 


2029 


8.48 


3.37 


8.03 


0.11 


0.07 


0.28 


0.68 


MC 47 


307 


363 


1.88 


0.64 


1.53 


0.12 


0.14 


0.47 


0.72 


MC 57 


337 


509 


2.25 


0.81 


1.93 


0.12 


0.13 


0.45 


0.74 


MC 64 


524 


826 


3.55 


1.65 


3.93 


0.090 


0.24 


0.66 


1.03 


MC 71 


594 


819 


3.83 


0.55 


1.31 


0.29 


0.11 


0.39 


0.34 


MC 90+91 


204 


382 


1.48 


0.87 


2.07 


0.071 


0.10 


0.37 


0.95 



Note. — Col.(l): Object name. Cols. (2) and (3): IRAS co-added flux densities at 60 and 100 fim, respectively. 
Col. (4): IR luminosities determined from cols. (2) and (3). The adopted distance is 51.8 kpc. Col. (5): Observed 
5 GHz flux densities. Col. (6): LC photon production rates estimated from col. (5) when the distance is adopted 
as 51.8 kpc. Col. (7): Ratios of IR luminosity to LC photon production rates. Col. (8): Eb-v estimated from the 
observed Balmer decrement. Col. (9): e determined from Figure 1. Col. (10): Estimated / from equation (2). 
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Table 5 

Summary of determined / for sample galaxies 



Galaxy 


/ 


12+log(0/H) 


V/Vmw 


(1)' 


(2) 


(3) 


(4) 


Galaxy 


0.45 ±0.21 


8.7 


1.0 


M31 


0.38 ±0.14 


9.0 


1.2 


M33 


0.41 ±0.37 


8.4 


0.6 


LMC 


0.74 ±0.22 


8.37 


0.2 



Note. — Col.(l): Galaxy name. Col. (2): Mean / 
with the sample standard deviation. Col. (3): Metal- 
licity from van den Bergh (2000). Col. (4): Dust- 
to-gas mass ratio normalized by the Galactic value, 
6 x 10~ 3 . 
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APPENDIX 

INFRARED EXCESS (iRE) 

The IR excess (IRE) is often examined for H II regions. We examine the relation between the LCE and IRE. Especially, 
we relate / to IRE directly here. A concerned discussion is found in Hirashita ct al. (2001). 
Mezger et al. (1974) have defined the IRE as 

IRE = — ^ , (Al) 

where h is the Plank constant and J^Lya denotes the frequency at the Lyman a line. From equation (Al), we obtain 
£ir,7/Al C 50 = 4.23 x 10 ~ 2 1 RE. Therefore, equation (2) is reduced to 

= 0-44 + 0.56. 

0.28 + 0.2AIRE ' v ' 

In Figure A5, we show the relation between / and IRE derived from equation (A2) for various e. 

The solid, dotted, and dashed lines correspond to t = 1.0 (in the Galaxy), 0.7 (in M31 and M33), and 0.4 (in LMC), 
respectively. In addition, we also present the dash-dotted line of e = for comparison, which corresponds to the case 
that no nonionizing photons are absorbed by dust. According to Spitzer (1978), the probability of producing Lyman a 
photons is two-thirds per every ionization-recombination process of hydrogen. Thus, when the observed IR luminosity of 
an H II region is explained by only the luminosity of Lyman a photons produced in the region (/ = 1 and e = 0), IRE 
will be 0.67. 

Antonopoulou & Pottasch (1987) have reported 4 < IRE < 40 for the Galactic compact H n regions. The mean IRE of 
their sample is about 12. For such high IRE regions, we expect that / is about 0.3 or less. When we estimate the SFR of 
such regions, therefore, we should not use the Ha or thermal radio luminosities but use the IR luminosity as an indicator 
of the SFR. 
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E-B-V 

Fig. 1. — Average extinction efficiency for nonionizing ultraviolet photons, e, as a function of color excess, Eb-v- The solid lines arc made 
by using the Galactic extinction curve (Seaton f979). The LMC extinction curve (Howarth f983) results in the dashed lines. The thick lines 
denote the results of cases that dust albedo is set to 0.5. The thin lines are the results of no scattering. 
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I I I I I I I I J I I I I I I I I I I 



.8 — 



.6 — 



.4 — 



.2 — 



LMC )< 







i i i i I i i i 



M33 



Milky Way 



M31 



I I I I I I I I I 



7.5 



8 



8.5 



9.5 



12+log(0/H) 



Fig. 2. — Fraction of Lyman continuum photons contributing to hydrogen ionization, /, as a function of metallicity. The data of 12+log(0/H) 
are taken from van den Bergh (2000). The error bars denote the sample standard deviation of the mean presented in Table 5. 
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-2-10 1 



log (D/D MW ) 

Fig. 3. — Fraction of Lyman continuum photons contributing to hydrogen ionization, /, as a function of dust-to-gas mass ratio, V. The 
dust-to-gas ratio is normalized by a typical Galactic value, T>mw = 6 x 10~ 3 (Spitzer 1978). The open square represents the case that T> of 
M31 is taken from Issa et al. (1990). The error bars denotes the sample standard deviation of the mean. Two solid lines represent predictions 
by the model that the dust optical depth over the ionized region is proportional to T> like equation (5). 
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log ( D/ D m ) 

Fig. 4. — Correction factor for star formation rate by the LCE effect as a function of dust-to-gas ratio. The correction factor is 1//. The 
below and above lines are the case of x = 1 and 2 in equation (5), respectively. The correction factor is expected in the shaded area between 
two curves. 
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Fig. A5. — Fraction of Lyman continuum photons contributing to hydrogen ionization, / as a function of infrared excess, IRE. The solid, 
dotted, dashed, and dash-dotted lines are calculated by equation (A2) with e = 1.0, 0.7, 0.4, and 0, respectively. 



